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Abstract 

We show that thermal effects play a relevant role in the determination of the 
electron-phonon couplings based on the intensity of the phonon-shakeup peaks 
in photoemission spectra. In particular, we re-consider the determination of 
the electron-phonon couplings of fullerene based on a fit of the photoemission 
spectrum of Cg in gas phase (1). We show that taking thermal effects into 
account reduces the obtained couplings by approximately 10%. 



INTRODUCTION 



Many of the interesting phenomena observed in electron-donor doped fullerides are 
related to electron-phonon coupling of the partly filled band derived from the lowest 
unoccupied molecular orbital (LUMO) of C 60 . These couplings, in turn, are mainly 
of molecular origin, involving the optical phonons derived from the molecular modes 
of A g and H g symmetry. Several different determinations of the electron-phonon cou- 
plings of the LUMO are available on the market, both from calculations (2-5), and 
derived from experimental data (1). All calculations based on density- functional the- 
ory (DFT) invariably turn out much smaller couplings than those obtained through 
fitting the experimental photoemission spectrum (PES) of gas-phase Cg (see for ex- 
ample Table GJ. This spectrum is affected by the electron-phonon shakeup associated 
to relaxation after the sudden ejection of the extra electron from the LUMO. That 
fit was realized on the basis of a zero-temperature model for the electron-phonon 
coupled system. The authors of Ref. 1 estimate the vibrational temperature of the 
sample ions of the order of 200 K. We claim that even such a low temperature could 
affect significantly the spectrum, in a way not unlike an increase of the couplings. 
We conclude that a determination of the couplings via a model including realistic 
thermal effects must turn out smaller couplings than a zero-temperature model. 



Mode® Exp. freq. (cm- 1 ) DFT Energy (meV) g Aj PES (1) g Aj DFT (5) 



AW " 496 ~~~ ~"~ 0.0 " 0.157 

A® 1470 0.851 0.340 

HW 271 32 0.824 0.412 
437 53 0.941 0.489 
710 89 0.421 0.350 
774 97 0.474 0.224 
1099 139 0.325 0.193 
Hf 1250 158 0.197 0.138 
Hp 1428 180 0.339 0.315 
1575 197 0.376 0.289 



Table 1: Experimental frequencies of the vibrational modes of fullerene linearly 
coupled to the t\ u LUMO. Computed energies and dimensionless coupling constants. 

2 THE MODEL 



The standard formulation of the linear electron-phonon model (1,6,7) describing the 
JT coupling of the orbitally degenerate t lu molecular triplet with the vibrations of 
symmetry A = A g , H g is the following: 

H = H + H vih + H c _ v , [2.1] 
H o = e^ a Lcma, [2.2] 

ma 

H vih = W^^{Pl + Q%) , [2-3] 

H c - V = k A ]T C^^Y^QKj^KjQm^naCm'a- [2-4] 
a fimm' j 

Here, Qj> = (b^ + bj^ / V% is the dimensionless normal-mode vibrational coordi- 
nate (in units of the natural length scale of the harmonic oscillator), and Pj^ the 
corresponding conjugate momentum for mode j of symmetry A (either A g : 2 modes, 
or Hg. 8 modes). & ma is a fermion operator creating a spin-a electron in orbital 
m = 0, ±1 of the LUMO. to, to' and /x label components within the degenerate 
multiplets, for example according to the C 5 character in the X D D 5 D C 5 group 
chain (8,9). In particular, /i — 0, ±1, ±2 for H g modes and /i — only for A 9 modes. 
C^, are Clebsch-Gordan coefficients (9) of the icosahedral group X that couple two 
ti„ tensor operators and a A tensor operator to a global scalar A operator (in practice 
these coefficients equal spherical Clebsch-Gordan coefficients for coupling two / = 1 
angular momenta to L = for A = A g and L = 2 for A = H g ). Numerical factors 
k A = \/3, k H = y3/2 (which could otherwise be re-absorbed into the definition of 
g\j) are included to make contact with previous notation (5,7). g\j are precisely the 
electron-phonon coupling parameters whose determination is addressed in this work. 
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Figure 1: The schematic "Frank-Condon" structure of the PES experiment, as seen 
from the point of view of any given normal-mode coordinate Q. Along each normal- 
mode direction, the displacement of the minimum of the ionized configuration with 
respect to that of the neutral configuration is proportional to the coupling of that 
mode. 

The Hg-mode coupling in the model |2.1j is of the Jahn- Teller (JT) type: exact 
solution are known only for very large or very small couplings. The couplings of 
fullerene appear to be in an intermediate regime, for which only numerical solutions 
are available. Indeed, the authors of Ref. 1, applied Lanczos diagonalization to the 
Hamiltonian \2.1\ on a truncated harmonic basis, to obtain its ground state. The 
analysis of this ground state in terms of the harmonic ladder of oscillators of the non- 
JT neutral C 60 provided the intensities to reconstruct the spectrum. At any finite 
temperature one should in principle repeat the analysis not only for the JT ground 
state, but also for a number of low-lying excitations and sum the resulting spectra, 
each multiplied by the corresponding Boltzmann weight. An attempt to follow this 
route at any finite temperature such that k^T is comparable to even the lowest 
quanta of vibration (« 30 meV) is beyond the possibilities of nowadays computers, 
due to the large number of oscillators, and the associated rapidly increasing density 
of coupled states: a reliable evaluation of the spectrum would require the accurate 
determination of the wavefunction of several hundred thousand of the lowest states 
in a Hilbert space of several million. 

We propose here a more efficient method, based on a related much simpler model, 
for which an analytical solution allows to carry out the calculation even to fairly high 



temperature, with mild computational requirements. The basic observation here is 
that in the experiment under examination the final states (neutral Cgo) involve simple 
harmonic ladders of phonons. These are excited by an amount which is related to the 
displacement of the adiabatic minimum of the initial C^ with respect to that of the 
final Cqq, as illustrated in Fig. ^ The intricate details of the initial JT state should 
play only a minor role in this spectrum, the main thing that needs to be made right 
being the displacement between the two minima. We replace therefore each of the 
H g modes (with the complication of the full JT Hamiltonian [2.1] ) with a set of five 
effective A g modes, with the same frequency as the original fivefold-degenerate mode, 
and a suitable coupling g^. It is important to introduce 5 degenerate oscillators per 
mode, to reproduce correctly the increasing density of oscillator states with energy, 
which is crucial to describe thermal effects. The advantage of A g nondegenerate oscil- 
lators is that their wavefunction expansion on a displaced basis is known analytically. 
Note that this approximation could not possibly work in a context (such as that of 
the PES from neutral C§q (10, 11)) where the final states are affected by JT. 

Briefly, starting from an initial state \i) = \v) with v quanta in the initial state, 
the required probability amplitude to end up into a final state |/) = T~ a \w) (T~ a is 
a space translation operator and the distance a = \ g&) with w quanta is 



</|aLl<> = MT a M = ^^"Cl ^ ) [2-5] 
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where L*(z) are the Laguerre polynomials of degree r: L^(z) = - r f £- F yz r+x e 
The spectral function is then computed according to Fermi's golden rule in the sudden 
approximation: 

h{hv) = -^{ftV^Mhv - {E f - Et)), [2.6] 

where Ei and Ef have trivial harmonic expressions. Assuming a Boltzmann popula- 
tion of the initial states 
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[2.7] 



the final spectrum is obtained through the following thermal average: 

I(hu)=Y,P(i)Ii(H- [2-8] 

i 

In the problem of C^q at hand, the total number of oscillators is 5 x 8 + 2 = 42. 
To combine the spectral contributions of individual oscillators, the collective matrix 
element in |2.6j is the product of the matrix elements of the individual oscillators, 
while the position (Ef — E^ of the peak is of course the sum of the energy change of 
the individual oscillators: 

(f\Vj\t) = (w 1 ,w 2 ,...\V I \v 1 ,v 2 ) = (w 1 \T^\v 1 )(w 2 \T a '\v 2 )... [2.9] 
E f -Ei = (E W1 + E W2 + ...) - (E V1 + E V2 + ...) = [2.10] 
= (E Wl — E Vl ) + (E W2 — E V2 ) + ... 
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Figure 2: The best fit value of function of gn, according to a T = 

calculation. Arrows point to the values of the total g H according to the DFT (5) and 
PES (1) estimates. 

The delta-function peaks are finally broadened into Gaussians of 20 meV HWHM to 
account for the final resolution of the spectrometer. For more details on this simple 
model, the reader is referred to Ref. 12. 

To determine the value of the effective couplings g^., we fit the T = spectrum 
of a single H g mode (computed by exact Lanczos diagonalization of |2.1j ) with that of 
5 degenerate A g oscillators of the same frequency. By repeating this fit for a range of 
couplings, we obtain the curve of Fig. El We have then verified that the temperature 
dependence of the spectrum based on the zero-temperature fit is in very good accord 
to that obtained by exact diagonalization (which, for this single-mode case, can be 
carried out up to moderately high temperature). 

For treating many modes, it is necessary to take into account the "collec- 
tive" nature of the JT interaction. Therefore, instead of converting each individ- 
ual guj through the relation of Fig. |2l we consider the "total coupling" (13, 14) 

~g H = 9 2 Hj) > an d convert g H into a , through Fig. El We then multiply each 
individual gnj by the same ratio /9h^ to obtain the effective g e j^ couplings for each 
mode. For the DFT coupling reported in Tabled the ratio I~9h = 1-045, while for 
the PES couplings Wa I~9h = 0.859. On the basis of these couplings, we compute the 
T = spectrum, and in both cases we obtain a good accord with the exact Lanczos 
calculation. Figure E2 illustrates this comparison for the PES couplings. 
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Figure 3: T = comparison between the exact JT Lanczos spectrum (solid line) and 
that obtained by the simplified effective- A modes model (dashed line). Comparison 
done for PES couplings (1), Gaussian broadening of 10 meV HWHM, including only 
H g modes. 

3 RESULTS AND DISCUSSION 

We apply now the effective- A modes model to the computation of the PES at finite 
temperature. Figure 0] reports the thermal spectrum for realistic values k^T = 20, 
40 meV, computed on the basis of the DFT parameters (5), rescaled as discussed 
above. Clearly, even though thermal effects go in the direction of increasing the 
transfer of weight to the phonon excitation and reducing the distance to the ex- 
perimental spectrum, the agreement remains poor, thus confirming that the DFT 
parameters are indeed too small. There is not a satisfactory explanation yet for this 
failure of DFT, which is particularly surprising in light of the accurate values for the 
couplings of the HOMO of C 60 obtained by the same method (11). 

We then verify thermal effects on the spectrum based on the couplings param- 
eters fitted on the PES spectrum using the T = model (1). Figure El reports the 
thermal spectrum for k^T = 0, 20, 40 meV. Clearly, for T = 0, the accord with exper- 
iment is good. As T is raised, more and more weight moves to the phonon satellites, 
making the accord less good. To restore a better accord at finite temperature, smaller 
couplings should be considered, reduced by perhaps 10%. For a precise evaluation of 
the amount of this reduction, a reliable estimate of the temperature of the sample is 
needed. Unfortunately, the intensity of the "pre-edge" anti-Stokes feature is only a 
rather indirect measure of temperature, since this intensity depends crucially also on 
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Figure 4: All-modes PES computed based on the effective- A modes model at fceX '=20 
and 40 meV, for the DFT couplings of Ref. 5. The poor accord with experiment 
indicates that, even accounting for thermal effects, the DFT couplings are indeed too 
small to account for the observed transfer of spectral weight to the phonon satellites. 

the value of the coupling of the involved low-frequency phonon mode. (This is clearly 
illustrated in the comparison of the k^T = 40 meV dashed curves of Figs. 0]and[5l). 
To obtain more reliable couplings, an independent estimate of the vibrational tem- 
perature of these ions would be necessary. It would then be possible to re-compute 
the electron-phonon parameters by means of a new fit, including thermal effects. 
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